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The effect of contact torques on porosity of cohesive powders 



o 
o 

(N 



o 



i 

c 

o 
o 



> 
o 



m 
o 

o 



o 
o 



13 



Guido Bartels 1 , Tamas Unger 1 ' 2 , Dirk Kadau 1 , Dietrich E. Wolf 1 and Janos Kertesz 2 



Abstract The porosity of uniaxially compacted cohesive 
powders depends on the applied stress (including gravity). 
The case, where these stresses are weak, is considered. 
The compaction results in a porosity which is a function 
of sliding, rolling and torsion friction. By contact dynam- 
ics simulations it is shown that the influences of contact 
torques (static rolling and torsion friction) on the porosity 
are significant and approximately additive. The relevance 
for nano-powder pressure sintering is discussed. 
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Introduction 

The behavior of granular packings under external load is 
governed by particle properties (e.g. roughness, elasticity, 
cohesion) as well as by the geometrical structure of the 
packing (e.g. connectivity, orientations of contacts). When 
compacting loose granular material, density and connec- 
tivity increase until a static state is reached, where the 
material withstands the external pressure. Such jammed 
states are currently widely studied 

For noncohesive materials the porosity in such a 
jammed state depends essentially on the deformability of 
the grains under the fixed external load. If the particles 
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were rigid, the porosity could not be reduced by upscaling 
the external load (including gravity) , but only by shaking 

HE]. 

However, for cohesive powders, static states with much 
higher porosity are possible. They are due to stabilization 
mechanisms, which involve the cohesion force as an in- 
trinsic force scale. Therefore, these states are not stable 
for arbitrary upscaling of the external load, but only up 
to a threshold depending on the porosity. Exceeding this 
threshold leads to further compaction of the powder. In or- 
der to avoid confusion we call these static states "blocked" 
instead of jammed. 

This paper considers stabilization mechanisms of 
blocked states in three dimensional cohesive powders of 
rigid particles. The relative motion of two solid spheres 
has six degrees of freedom, three translational charac- 
terized by a velocity vector v with one normal compo- 
nent and two tangential components (sliding modes), and 
three rotational characterized by an angular velocity vec- 
tor uj, again with one normal component (torsion mode) 
and two tangential components (rolling modes). If any of 
these modes are damped or blocked they represent specific 
dissipation or stabilization mechanisms, respectively. Well 
known examples are the static and sliding friction: the 
former stabilizes a contact against small tangential forces 
(blocking), the latter dissipates kinetic energy (damping) 
in case of sliding. 

All modes of relative motion can be blocked: In addi- 
tion to static friction also resistance against separation, 
rolling and torsion must be taken into account, i.e. the 
contact can exert also cohesion force, normal and tangen- 
tial torques (i.e. torsion and rolling torques, respectively) 
in order to inhibit relative motion. To what extent these 
particle interactions stabilize pores in a cohesive powder 
is the subject of this paper. 

Extending previous work in two dimensions we 
introduce torsion friction in addition to rolling friction 
and cohesion into contact dynamics simulations of a uni- 
axial compression process. Of course, two perfectly rigid 
spheres, if they existed, would only have a contact point, 
which could neither exert a rolling nor a torsion torque 
on the particles in contact. Even the Coulomb- Amonton- 
daVinci friction law would not be justified for such an ide- 
alized point contact. By contrast, we consider rigid spheres 
here only as a geometrical idealization of real particles 
and do not take their contact areas into account explic- 
itly. Implicitly, however, the finite size of the contacts is 
responsible for the various kinds of friction we consider. 
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It is believed that torsion and rolling resistance are 
of little importance in noncohesive granular assemblies of 
spheres, where indeed in statics the contact-torques vanish 
jjj or are very weak [H]. Our work was motivated rather 
by powders with grain size smaller than lOOnm. The field 
of nano-powders [5j attracts a lot of scientific and in- 
dustrial interest due to their different material properties 
compared to assemblies of larger grains or bulk materi- 
als. Many concepts developed for regular granular media 
can be applied also to nano-powders, but additional as- 
pects have to be taken into account as well, such as strong 
cohesion (due to van der Waals forces) and sinter-neck 
formationjlOj between the grains which make the ques- 
tion of blocking torsion and rolling modes relevant. Little 
is known so far about the way in which a sinter-neck re- 
sists rolling or torsion. It is plausible, however, that on a 
sufficiently short time scale (where creep can be neglected) 
torque thresholds have to be exceeded to break a sinter- 
neck and induce relative motion between the particles. 
Therefore, lacking well proved contact laws on the nano- 
scale, we assume the simplest kind of threshold dynamics 
for sliding, rolling and torsion friction. 
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Model 

We consider rigid spherical particles of identical sizes, but, 
as mentioned above, force and also torque transmission is 
allowed at contacts. The total normal force F n between 
two grains has two constituents: first the attractive part, 
which is a constant cohesion force F co h and second the 
force -Fexc due to the excluded volume constraint: 

Fn ^cxc -^coh ; (1) 

where the repulsion (attraction) is denoted by positive 
(negative) sign. Depending on external forces acting on 
the particles the constraint force -F e xc can take any positive 
value. As a consequence the total normal force can be 
repulsive F n > (arbitrarily strong) or attractive > 
F n > — i^ co h as well and is in this latter case limited by 
the cohesion force. Thus F n inhibits detachment as long 
as the pulling is weaker than the cohesion force Fooii- If 
-Fcoh is exceeded, the contact starts to open, but breaks 
only when the work done by the pulling force exceeds the 
cohesion energy (F con times cohesion range). More details 
of this implementation of cohesion can be found in JI] . 

The tangential force, as well as the normal and tangen- 
tial torques are responsible for blocking the sliding, torsion 
and rolling modes, respectively. Their implementation is 
described next. 

For the tangential contact force the Coulomb friction 
law has been modified to include the influence of cohesion: 
If the relative velocity of the two surfaces is zero (sticking 
contact) the magnitude of the friction force (|Jt|) can be 
any value up to the threshold ft, max, while in the case 
of sliding |F t | = ft, max, and its direction is opposite to 
the relative velocity (Fig. The maximal friction force 
is proportional to the normal force including F coa : 

-Ft, max = n(F n + F con ) , (2) 

thus the threshold value vanishes when the contact opens 
(Fn = — fcon)- 
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Fig. 1. The graph represents Coulomb's law of friction: A 
sticking contact (with relative tangential velocity vt = 0) can 
bear any tangential force with absolute value up to ft, max- If 
sliding occurs (v t ^ 0) the magnitude of the tangential force is 
ft, max while its direction is opposite to the sliding velocity. 
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Fig. 2. Contact torques applied to suppress rolling and tor- 
sion. The rules are similar to those of the tangential force. 
The rolling torque (torsion torque) is coupled to the tangential 
(normal) relative angular velocity oj t . rc i (w n , re i). For definition 
of the threshold values T„, max and Thmax see the text. 



The contact laws for the normal and tangential torques 
are chosen in analogy to the Coulomb friction law (see Fig. 
|2J. The threshold values for the torques T n max and Tt imax 
are defined as force times length, where the force-scale is 
again given by (F n +F con ) and the radius of the sinter- neck 
determines the relevant length-scale (Fig. 0) . This length 
makes sense only if it is below the particle size, and in our 
study we choose it within the range between zero and the 
effective radius at the contact: r c s = (1/Ri + I/-R2) , 
which for two identical spheres is half the particle radius. 
The maximal normal and tangential torques in the model 
are thus given by: 



-Fn,max Mn(fn ~\~ fcoh)^eff 
.max — 



(3) 
(4) 



where the dimensionless parameters /i n and fit make it 
possible to control the strength of blocking, similar to the 
friction coefficient /z in Eq. 

This model is certainly oversimplified, and the physics 
of nano-particles might be better described by more so- 
phisticated, yet unknown microscopic models. Moreover, 
in our model rolling, torsion and sliding modes are inde- 
pendent of each other. In general this is not the case as 
demonstrated for an ordinary disk on a flat surface [511121 
113) . where the coupling between torsion and sliding fric- 
tion reduces both of them. In this case assuming indepen- 
dence would overestimate friction. If this is generally true, 
the porosities obtained in our simulations are expected to 
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Fig. 3. Schematic figure of a sinter neck between two particles. 

be upper bounds. With some caution we can therefore use 
our simplified model to analyze the effect of rolling and 
torsion friction on the porosity. 



3 

Computer Simulations 

The discrete element simulations presented here were per- 
formed with the contact dynamics method [14115116) . 
where the particles are modeled as non-deformable spheres 
and the interaction between particles is determined by 
constraint conditions. The application of constraint forces 
and torques and an implicit time stepping scheme make 
this method especially suitable for the implementation of 
threshold dynamics and for the proper treatment of static 
forces and torques in blocked or jammed states. 

More about our algorithm can be found in [1711181119) . 
These works contain the description of the basic 3D algo- 
rithm for ordinary granular materials, an analysis of the 
method and its extension to the case of nano-particles. 

The system under consideration consists of spherical 
particles with the same radius R. Starting from a low den- 
sity state the compactibility is measured for various values 
of /it and // n , while the coefficient of friction, \i = 0.3, is 
kept constant. 

The initial configuration is prepared by ballistic de- 
position 20 2"T|lr)] : Particles fall vertically towards a hor- 
izontal plane, one by one with randomly chosen x — y- 
coordinates. As soon as a falling particle comes closer 
than a capture radius r capt to the deposit or substrate, 
the contact is established immediately, and the particle 
sticks irreversibly. Then the next particle is dropped. Here 
we choose r capt = 3R. As a result one obtains deposits of 
very low density with tree- like structures. Using this type 
of initial configuration is motivated by filter processes ex- 
tracting nano-particles out of a gas flow [221123) . 

We simulate the uniaxial compression of this ballistic 
deposit by a piston moving along the z-direction towards 
the bottom plane, with periodic boundary conditions in 
x- and y-direction. Gravity is neglected. The compression 
is caused by a constant pressure F/L x L y on the piston 
(Fig. [3}. The simulation ends, when the system reaches 
a blocked state, i.e. when the piston comes to rest. The 
system has size L x = L y = 25i?, contains 1015 parti- 
cles and the mass of the piston is chosen to lOOOpi? 3 , 
where p denotes the mass density of the particles. The 
pressure on the piston is F co h/(400i? 2 ) in all cases. This 
is a very weak compacting pressure, as the typical dis- 




Fig. 4. a) Initial arrangement with 1015 ballistically deposited 
particles, b) Final configuration of the system compressed by 
the external force F. 

tance between branches of the ballistic deposit is of the 
order of 5R, so that the typical load on a single branch, 
i r ' co h25/400 = F co h/16, is much smaller than the intrinsic 
force scale given by the cohesion force .F co h. 

The mesh of values (/it, A* n ), for which we determined 
the final porosity, is indicated in FigEJ Assuming a smooth 
dependence of the porosity on the friction parameters, we 
can estimate the error bars from Fig0 without the need 
to do more than one run for each data point. 
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Results 

Based on the final position of the piston we measured the 
porosity E of the blocked states: 



i.e. the relative free-volume in the system (here V de- 
notes the volume between the bottom and the piston). 
Fig. [5] shows that contact torques have significant effect 
on the compactification: Whereas without rolling and tor- 
sion friction the final porosity is Eq rj 54% for the weak 
compaction pressure we considered, porosities as high as 
82% are stable for /z t = = 1- 

The porosity added due to rolling and torsion friction, 

Ei(pt, /x n ) = E - E , (6) 

saturates in the region where the coefficients /i t and /z n 
are larger than 0.4. More than Ei(l, 1) ss 28% cannot be 
achieved based on the contact torques, not even if one 
increases the threshold values far beyond the physically 
meaningful range. The maximum porosity is less than the 
one of the initial configuration: Suppression of rolling and 
torsion degrees of freedom alone does not suffice to avoid 
compaction completely. 

Qualitatively, rolling friction alone (/i n = 0) as well as 
torsion friction alone (fit — 0) have similar effects on the 
porosity. Quantitatively we find, that -Ei(l,0) » 18% is 
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Fig. 5. Final porosity of the compacted system obtained for 
various values of rolling (pit) and torsion resistance (/x n ). 
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Table 1. The influence of the normal and tangential contact 
torques is shown on the stress transmission along the direc- 
tion of the uniaxial compression. af z and a^ z are containing 
only contributions of the compressive and the tensile contacts 
respectively. 



about twice as big as £i(0, 1) « 8%. The reason is not 
clear, but it is intriguing to notice that rolling friction /x t 
suppresses two degrees of freedom, while (j, n suppresses 
only one. The ratio E\{fi, 0)/Ei(0, p) » 2 is roughly inde- 
pendent of jl. 

An interesting property of the function fii(/it, /i n ) is, 
that it can be very well represented by the following sum: 

£i(/i t ,/v) =£i(/it,0) + £i(0,Mn) • (7) 

This shows that rolling and torsion friction contribute in- 
dependently to the porosity. The difference \AE\ between 
the two sides in Eq. (J7J) is less than about 3%, and AE 
fluctuates around zero with no apparent systematics. 

Finally we would like to discuss, how the static re- 
sponse of the system to the external load differs with and 
without contact torques. In the absence of contact torques, 
strong tensile contact forces appear and seem to be crucial 
for stabilizing the compressive force lines against buckling 
in the porous system. This pore stabilization mechanism 
has also been found in two-dimensional systems 6 . Of 
course, the contribution of the compressive forces to the 
macroscopic stress tensor overcompensates the one of the 
tensile forces in order to balance the external load (table 

m 

In the case where contact torques are allowed, rolling 
and torsion friction already stabilize the force lines against 
buckling, before significant tensile forces develop. This 
does not mean that cohesion is superfluous in this case 



as it is also responsible for the enhanced threshold values 
of the contact torques (Eq. J3J|). 

We measured the macroscopic stress tensor ( |24| ) : 

fX'-:^ fX':'; % • (8) 

c<+) c(-) 

separately taking only the contacts under compression 
(crtj) or the ones under tension (cH) into account. (Here 
lj is the jth component of the interstice vector connect- 
ing the centers of the particles at contact c.) The zz- 
components of the stress tensor are shown in the table , 
where the resulting stress a zz = F/L x L y is of course de- 
termined by the pressure on the piston, which is the same 
for all cases and is taken as unit stress. In the zero-torque 
case more than twice of the external pressure is provided 
by the contacts under compression, because a sufficient 
amount of tension must be allowed in the system in order 
to stabilize the pores. This internal counter-stress is di- 
minished significantly by rolling friction, whereas torsion 
friction alone has little effect on the partial stresses <r+, 
and a~ z . 
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Conclusion 

We presented 3D contact dynamics simulations of cohe- 
sive powders and studied the porosity under weak uniaxial 
compression. We focused on the effect of contact torques 
which suppress relative torsion and rolling of the adja- 
cent particles. It was found that the presence of contact 
torques has enormous impact on reducing the final den- 
sity of the system and in addition the characteristics of the 
stress transmission is altered: Without torsion and rolling 
friction strong tensile forces develop and play important 
role in the mechanical stabilization. These tensile contacts, 
however, cancel a large part (more than half) of the pres- 
sure exerted by compressive forces and only the remaining 
part is utilized to resolve the external load. This situation 
is changed by switching on contact torques, which reduce 
tensile forces significantly in the system. In that case the 
stress due to the compressive contacts corresponds ap- 
proximately to the external load and only a small part of 
this stress is "wasted" to overcompensate tensile forces. 

A remarkable feature is the additivity found in the 
porosity, i.e. the porosity is well represented as the sum of 
independent contributions of the torsion and rolling fric- 
tion. Our results suggest that each relative motion mode, 
if it is suppressed, results in an additional free volume in 
the system independent of the other modes. Whether this 
picture holds also for translational degrees of freedom (i.e. 
altering the Coulomb friction coefficient or the cohesion 
force) is a subject of further investigation. 
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